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A unified viscoplastic constitutive relation based on crystallographic 
slip theory is developed for the deformation analysis of nickel-base face- 
centered cubic superalloy single crystals at elevated temperature. The 
single crystal theory is then embedded in a self-consistent method to de- 
rive a constitutive relation for a directionally solidified material comprised 
of a polycrystalline aggregate of columnar cylindrical grains. One of the 
crystallographic axes of the cylindrical crystals points in the columnar 
direction whilst the remaining crystallographic axes are oriented at ran- 
dom in the basal plane perpendicular to the columnar direction. These 
constitutive formulations are currently being coded in FORTRAN for use 
in nonlinear finite element and boundary element programs. An experi- 
mental program to determine the biaxial tension-torsion behavior of PWA 
1480 single crystal tubular specimens at 1600°F is also underway. 


INTRODUCTION 

This paper represents a first quarterly progress report on a program 
to develop anisotropic constitutive equations for use in modeling the mul- 
tiaxial viscoplastic stress-strain response of single crystal and directionally 
solidified gas turbine alloys at elevated temperature. Two approaches are 
being pursued. The first approach consists of modeling the anisotropic 
response from a macroscopic continuum point of view, whilst the second 
approach consists of taking a micromechanics viewpoint using crystal plas- 
ticity concepts. Tubular specimens of the single crystal superalloy PWA 
1480 are currently being machined. Tension-torsion experiments on tubu- 
lar specimens of PWA 1480 will be conducted at 1600°F to provide a data 
base for exercising the theoretical formulations. The anisotropic visco- 
plastic theories will be incorporated into a nonlinear finite element code 
since the non-uniform stress distribution in the tubular specimens will 
require the solution of a boundary value problem for data reduction pur- 
poses. 


65 



SINGLE CRYSTAL ANALYSIS 

An analysis of single crystal superalloys undergoing steady state 
creep deformation was presented by Paslay, Wells and Leverant [l] in 
1970 using a theoretical formulation based on crystallographic slip theory 
of face-centered cubic materials. In 1971 the theory was applied by 
Paslay, Wells, Leverant and Burck [2] to describe the creep behavior of 
single crystal nickel-base superalloy tubes under biaxial tension. Steady 
state creep formulations suitable for the analysis of single crystals were 
used by Brown [3] in 1970 and by Hutchinson [ 4 ] in 1976 to predict the 
behavior of polycrystalline materials whose aggregate consists of randomly 
oriented single crystal grains. Recently, Weng [5] has developed a single 
crystal creep formulation which accounts for transient (primary) as *veli 
as steady state (secondary) creep. However, in order to describe the 
combined plastic and creep behavior of polycrystalline materials, Weng 
combines the rate- independent and rate-dependent components of crystal 
behavior in such a way that each component is governed by a separate 
constitutive relation. The averaging of the single crystal creep relations 
to obtain the overall macroscopic creep response of the polycrystalline 
aggregate is easily accomplished by using -Kroner's self-consistent method 
[6J. In a general analysis the constitutive relations for the overall mac- 
roscopic plastic response of the polycrvstalline aggregate must be obtained 
using Hill's self-consistent method [7]. Kroner's method for calculating 
the macroscopic creep properties of a polycrystalline aggregate of single 
crystals is explicit in nature; but Hill's method for calculating the macro- 
scopic plastic properties is implicit in nature and requires lengthy itera- 
tive computations. 

In the decade of the seventies the creep and plastic responses of 
materials were combined into unified viscoplastic formulations [8J. These 
formulations differ from steady state creep theories by introducing history 
dependent state variables to account for primary creep and plasticity. 
Most of these unified theories exhibit an elastic response under instan- 
taneous deformation. This instantaneous elastic response occurs 
because the inelastic strain rate is assumed to depend only on stress, 
state variables and temperature and not on the rates of these variables. 
These unified theories may, therefore, be integrated in time by means 
of an explicit Euler forward difference method. Macroscopic properties 
of a polycrystalline aggregate, comprised of single crystal grains which 
are assumed to deform according .to a unified viscoplastic slip process, 
can therefore be obtained by means of Kroner's explicit self-consistent 
method. 

In the unit cell of the face-centered cubic crystal shown in Figure 1 
we denote byjJ;a unit vector in the i** slip direction (say of type <110> ), 
whilst xU is a unit vector in the normal direction to the slip plane (of 
type {111} ) of which constitutes a slip direction. The four octahe- 
dral {111} planes and the twelve corresponding <110^ slip directions 
(three on each plane) are shown in Figure L 

From the geometry of the unit cubic cell in Figure 1 the unit 
vectors are given by 
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2r=(i-i)/ fz , S*'(j +*)/& > »7=(Uj)/j2, Ii=(-itk)/^2 , 

* (t»j$)/t|2 , m, 0 »Cj*k)/j2 , «» = (t-j J/Ji , m,i= (Ul<)/Jj j 

with 

nv.j^x nj*(Uj + ]c)/j3 } n** (i+i-k)/^ , 

«7=st- C-i-i+iO/Js » 

where i; i» i are xmil vectors along the x,y,z crystallographic axes. 

Figure 2 shows a single crystal whose global axes are denoted by 

x , y , z and whose crystallographic axes are denoted by x,y,z. If Q 
denotes the orthogonal tensor which rotates the crystallographic (unstar- 
red) axes into the global (starred) axes, viz., x* * Q-,jX: , then the 

stress tensor jr and the strain rate tensor i in the crystallographic axes 
may be obtained from the stress tensor o - * and the strain rate tensor j|* 
in the global system from the usual transformation relations, 

JJ «<$.£*. Ct T and e«Q.fc*Q T . 

** ****** ******* (1) 


The assumption is now made that any of the unified viscoplastic 
theories discussed in Reference [8], when specialized to the case of shear 
deformation, is a valid constitutive relation in each of the twelve crystal- 
lographic slip directions. In the r** slip direction the resolved | shear stress, 
*t£ n » »s obtained from the relation 1 


= *!■.£•£»• for 


( 2 ) 


It is further assumed, in a manner analogous to the unified isotropic vis- 
copiastic models, that the applicable relation governing the inelastic shear 
strain rate in the r* slip direction is 

i r « K? {(iC-0|tC- ucr 1 ♦ 

+ ^ i (3) 

where K r and (with p and q = m,n,z) denote the drag stress and 

equilibrium (rest or back) stress in the r*' slip direction. The expression 
for Ytfwt is defined in equation (2). Terms such as denote, by an- 

alogy, expressions of the type 


TT 


r 

mz = 


. (T. Ir * 


(4) 


where Jgr denotes a unit vector perpendicular to the unit vectors m r and 
n r . The vector z r is in the slip plane containing the vector jinr and 
the vectors m r , n r , £** form an orthogonal triad for the slip sys- 

tem. In equation (3) the tensor otp^ represents the effect of the non- 
Schmid factors [9] upon the inelastic strain rate in the r* slip direction. 


6 7 



For example, the term containing represents the effect of the resolved 
stress, normal to the slip plane containing the r*' slip direction, on the 
inelastic strain rate in the r** slip direction* Such terms can represent 
the effect of a pressure dependent inelastic strain rate* The dominant 
term in equation (3) is the Schmid type term containing the expression 
iCin 5 estimates of the magnitude of the non-Schmid type terms con- 
taining the tensor have been given by Asaro and Rice in Reference 
[10]. * 

To complete the constitutive formulation it is necessary to specify 
the growth relations for the equilibrium and drag stress state variables* 
The equilibrium stress in the r* slip system may be assumed to evolve 
according to the evolution equation 

&k * ~ fpV^rl^ “ jJlWwl ** * (5) 

where and m are temperature dependent material constants. 

A simplifying assumption is to take 

4* - 4*- ^ f" p-'V"" 1 "' 1 ’ <«> 

K n r,,-*- . 

where the cumulative inelastic strain is defined by 




Equation (6) states that the ratio of the saturated equilibrium (back) 
stress, Le. the equilibrium stress ^f or continued inelastic straining 
under fast straining conditions ( Y r -» *o ), to the resolved shear stress in 
the same direction, is equal for all the non-Schmid systems. Moreover, 
the ratio is equal to the ratio of the saturated equilibrium stress to^ the 
resolved shear stress for the Schmid-type components, • 

Under continued inelastic straining the equilibrium stress saturates and 
Under fast straining conditions ( YJ.-* » ) equation (5) shows 
that the saturated equilibrium stress is given by the relation 

“n - 4- > (s) 

Tr-**° a 

since the thermal recovery term containing the material constant 

can be neglected for 9 r -+ oo . The material constants ^ must therefore 

satisfy the relation 

in _ Isa . * • (9) 

ft — f" ~ 

Assuming the constants ^ and and the limiting saturated value of 
the resolved shear stress are known for the Schmid-type terms, 

equation (9) determines the ratio of to fya for the non-Schmid 

terms. If it is further assumed that the initial hardening rate for the 
non-Schmid equilibrium stresses are all equal to the Schmid equilibrium 
stress hardening rate, then fort*l,2 and p^« «n>n,Z. . 
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The drag stress may be assumed to grow according to the evolution 
equation 

Kr S E > ( 10 ) 

Iwi 

in which the hardening moduli are given by 

. (ii) 

The hardening moduli defined in equation (10) account for the latent 
hardening effects observed in single crystal materials, and equation (11) 
is similar in form to that proposed by Hutchinson [11], Asaro [12] and 
Peirce, Asaro and Keedleman (13). Numerous forms for the hardening 

moduli hrk have been proposed in the literature and a review of harden- 
ing moduli may be found by consulting Havner's papers (cf. Reference 
[14]). The drag stress evolution equation does not contain thermal recov- 
ery terms, but these terms can easily be included in the formulation. 

The shear slip strain rates may now be resolved into the crystal- 
lographic system and summed for each slip system to obtain the inelastic 
strain rate tensor ^rith respect to the crystallographic axes in the form 

£ * + • ( 12 ) 

Finally, the stress rate tensor with respect to the crystallographic axes 
is determined from the relation 

z-gMi-i) > da) 

where £ is the anisotropic elasticity tensor for the face-centered cubic 
crystal referred to the crystallographic axes. 

The variables can now be updated in the Euler forward difference 

form: 

V-r.+Ut, g».£z.& , s*-s T .t.g , ( i4) 

where Afc is the current time increment. The process may then be re- 
peated by integrating equations (1), (2), (3), (4), (5), (10), (11.) for 
each time increment. 

The preceding discussion has focussed on slip which occurs on the 
t 1 11) octahedral planes in the <110> type directions of face-centered 
cubic nickel-base superalloys. Paslay, Wells, Leverant and Burck [1,2] 
also found that slip occurs under creep conditions on the {,111} planes 
in the <112) type directions. Cube slip on the crystallographic faces was 
also found to occur. For each different slip system the foregoing theory 
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is altered only by virtue of having different slip and normal vectors defin- 
ing the triad m, n and £ In general, the total inelastic strain rate may 
be written as the~sum, 

£ * H 

L»l 


where £t is due to {111} OlO^ type slip, £* is due to 4112^ 

type slip, c t is due to cube slip, etc. Such a combination was stated to 
be required to model primary creep behavior in Reference [l] by Paslay, 
Wells and Leverant. It is possible that {111} <112> type primary creep 
may evolve into {111} <110> type secondary creep as described by Lever- 
ant, Kear and Oblak in Reference [15]. In this instance it may be nec- 
essary to modify the theory so that the Xi, evolve with deformation, in- 
elastic strain rate and temperature according to evolution growth equa- 
tions similar to that employed for the equilibrium and drag stress state 
variables. This mixing of different slip systems to model the anisotropy 
of nickel-base superalloys was also stated by Ezz, Pope and Paidar [16] 
to be necessary in order to model tension-compression flow stress asym- 
metry observed in nickel-base superalloys. 


DIRECTIONALLY SOLIDIFIED ANALYSIS 


A model for directionally solidified alloys can be constructed by 
making use of a suitable self-consistent method to average the results of 
the single crystal viscoplastic constitutive theory. The directionally solidi- 
fied material consists of aligned columnar single crystal grains which are 
oriented at random in the basal plane perpendicular to the cylindrical 
growth direction. This random orientation of the grains produces a mat- 
erial with transversely isotropic properties. 

The directionally solidified material comprised of an aggregate of 
single crystal columnar grains may therefore be modeled in the following 
manner. We first choose a particular single crystal columnar grain and 
replace the aggregate of single crystal grains surrounding the chosen 
grain by a transversely isotropic material. The properties of this sur- 
rounding transversely isotropic material are found by averaging the prop- 
erties of the chosen single crystal grain about its cylindrical growth axis. 
It is then possible to relate the stress and strain increments in the single 
crystal gTain to those in the surrounding transversely isotropic matrix 
by means of the method proposed by Eshelby [17] in 1957. 

Viscoplastic formulations which exhibit an instantaneous elastic 
response can be integrated by means of an Euler forward difference 
method. In physical terms this integration process consists of letting 
the material creep at constant stress at the level jJ appropriate to the 
beginning of the increment for a time interval At . After the creep 
increment is completed the material is subjected to an instantaneous 
strain increment At - Ac , where A£ is the total strain increment during 
the time increment At and Afc is "the completed creep strain increment. 
The instantaneous application of the strain increment At-A£ induces an 
elastic stress increment given by Hooke's law in the form Ajr = JD:(A£-Ac) 
where £ is the elasticity tensor for the material. 
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Let 2 an ^ B denote the fourth rank cubic and transversely iso- 
tropic elasticity tensors for the single crystal and the surrounding matrix, 
respectively. In the cylindrical single crystal grain the inelastic strain 

increment is denoted by Ac and the corresponding quantity in the matrix 
is denoted by<Ac>. The quantity <Ac^ i s obtained from Ac by averag- 
ing Ac over ail angular orientations (viz. from 0 to 2rr ) in the basal 
plane perpendicular to the cylindrical growth axis. In the single crystal 
grain the constitutive relation has the form 

A<ry * Dijkl {tek\ “ A c kl ) 3 

whilst in the directionally solidified matrix the constitutive form is 
(Affjj) s Dyid ((AA K i) . (17) 

Given a known strain increment in the directionally solidified 

matrix, the object is to determine the corresponding inelastic strain 
increment Ac^ in the single crystal cylinder and then average this 
quantity by means of the relation 

<*C U > * to ** > (18) 

where 0 is the angle between the crystallographic axes x,y in the speci- 
men and the global axes x*,y* in the matrix, with the z,z* axes aligned 
in the cylindrical crystal's growth direction. The stress increment in the 
directionally solidified matrix is then determined by equation (17). 

The first step consists of determining the state of stress and strain 
in the cylindrical crystal grain when the matrix and crystal undergo creep 
for a time increment At . A preliminary step in this analysis consists 
of replacing the single crystal grain by a fictitious material which has 
elastic and inelastic properties which differ from those of the single 
crystal but in which the elastic properties are the same as that of the 
transversely isotropic matrix with elasticity tensor Jgf'. Eshelby's cutting, 
straining and welding operations [17] are now applied to the fictitious 
crystal grain. 

The fictitious cylindrical grain is now cut out of the matrix and the 
instantaneous shapes of the grain and resulting cylindrical hole in the 
matrix are maintained by appropriate equal and opposite surface tract- 
ions applied to the respective cylindrical surfaces of the grain and hole. 
The stress in the matrix is denoted by <£^ and that in the fictitious 
grain by £ , where £ is the stress state in the actual single crystal 
grain. From Eshelby's results, if the strain history in the matrix is 
homogeneous, the resulting stress in the actual and fictitious cylindrical 
grains will be constant throughout the cylindrical volume. This will be 
demonstrated subsequently. The fictitious cylindrical grain and the sur- 
rounding matrix are now assumed to undergo creep responses for a time 
increment At at their respective stress levels of £ and <£) . In the grain 
the resulting creep strain increment is denoted by Ae T and the corres- 
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ponding creep strain increment in the matrix is <Ae^ . If incremental 
surface tractions are instantaneously applied to the cylindrical surface 
of the fictitious grain so that it is elastically strained by an amount 
-a£ , it will regain the original size and shape which it had upon re- 
moval from the matrix prior to the creep response. If it is subsequently 
elastically strained by an amount (A£) through the application of a 
further set of incremental tractions applied instantaneously to the cylin- 
drical surface of the grain, it will fit back into the matrix from which 
it was removed. The fictitious grain and matrix now fit compatibly 
together, and the strain increment experienced by the grain is- where 


i|»l| * ^*1 ~ < iC H> ‘ 


However, a layer of surface traction exists on the cylindrical surface 
boundary between the grain and the matrix. This layer of surface trac- 
tion is given in magnitude by the relation ^ 

r\j * nj D^i AeJ, ) * *" n j Dijkl A£id j 

and can be removed by the application of an equal and opposite layer of 
surface body force of magnitude 

Af i * nj D^jkl Ajikl . (20) 

The displacement increment M; induced in the fictitious grain due to 
the application of surface tractions Af L over its cylindrical surface is 
given by 

a*! - 5S t ^ jO (21) 

where jr — jr # ) is the elastic Green’s function for the transversely 

isotropic matrix and the fictitious grain. Application of Gauss’ diver- 

gence theorem to the surface integral produces the result 

Au. C p * 5$ niDjn 

<22> 

The resulting strain increment in the fictitious grain due to this 

annihilation of the incremental surface traction built up during the creep 
response of the fictitious grain and the matrix is 

A£l f * 2 \ ’ J * * ( 23 ) 

where the Eshelby tensor 5, is defined by the relation 

s* - - i v { } . (24) 

The strain increment induced in the fictitious grain ir. order to make 
it fit compatibly in the matrix is given by " A(4 v j , whilst the additional 
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strain induced in the grain by removing the unwanted surface traction 
increment At;, is given by . At the end of the creep increment the 
total strain increment induced in the fictitious grain is given bvAz*-bfl- 
o. from equation (23), by SjMA*, . * The * “ g 'cre’Arenf^n 

the grain is therefore r J 

Affy • Iw )A» mn * (25) 

where I is the fourth rank identity tensor. Note that the strain incre- 
ment responsible for changing the size of the cylindrical grain from the 
size it had when it was placed back into the matrix is Aef:* S;:ki Mid . 
since the change in size is due to the annihilation of the surface trkction 
increment At; by the equal and opposite traction increment Af ; . 

Now consider the actual situation where the cylindrical grain has 
its own anisotropic elastic constant tensor with cubic symmetry, D c 
We remove the cylinder from the matrix, as before, and let the cylin- 
drical grain undergo a creep strain increment A£ in time At and the 
matrix undergo a creep strain increment <Ac> . If we now apply surface 
traction increments instantaneously to the cylindrical boundary so that 
the cylinder is elastically strained by an amount (At) - Ac it will 
fit back into the matrix from which it was removed. ^Moreover, since the 
matrix creep strain increment <Ac> i s the same as that in the problem 
With the fictitious cylinder, the cylindrical grain will have the same size 
and shape as the fictitious cylinder had when it was put buck into the 
^e . ar " n increment responsible for changing the of the 

cylindrical grain from. the size it had when it was placed back into the 
matrix to the size in its final configuration is Sirfmn A/&„,i . Hence if 

the actual cylindrical grain is strained by the increment S*U„ A/C 

.j r, Z have the same final size as the previously 

considered fictitious cylindrical grain with transversely isotropic elastic 
properties. The stress increment in the actual cylinder due to the strain 
increment + <Ac k <> - Ac k r is 

Ad\j = Djjw +<Ac kl }-Ac k | "J- . (26) 

If the stress increments in equations (25) and (26) are equal, the actual 
cylindrical gram which has the elasticity tensor D c appropriate to cubic 
symmetry and which undergoes an increment of transformation strain M 
may be replaced with the fictitious cylindrical grain with elasticity tensor 
2 equal to that of the matrix without upsetting continuity of dis- 
pl ace merits and tractions across the cylinder- matrix interface. Equating 
(25) and (26) shows that 

A^ij * * (27) 

If this, transformation strain increment occurs in a cylinder with elasticity 
tensor £ , the stress increment in the cylinder is equal to that which 

actually occurs in the single crystal cylinder with elasticity tensor D c 
undergoing a transformation strain increment A£-<Ac>. Substitution”of 
the expression for Ag, into equation (25) [or (26)] gives the stress incre- 
ment m the cylindrical single crystal grain at the end of the creep re- 
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( 28 ) 


sponse in the form 

AOij * 1* ( Ski**" ^PV* * ^ w ' r5 ^ 

* &e*U { ^ C tu " ^ -^ C tuV \ * 

According to the Euler foward difference integration method the 
crystal grain and matrix are now to be instantaneously loaded with the 
elastic strain increment — . To this end we first consider 

the single crystal to be replaced by a fictitious cylinder having the same 
transversely isotropic properties as the surrounding matrix. If the fic- 
titious cylinder undergoes a stress increment free uniform strain incre- 
ment At r , the final strain increment inside the cylinder is ( Sijkl — 

) heJi . A uniform strain increment - ^Ac> ma y n ? w ^ a P“ 

plied to the matrix and fictitious cylinder to produce the final strain 
increment of ( SijKl - Iijkl ) AtJ + <*£•)> - <Ac { j> . The resulting 

stress increment in the fictitious cylinder is 

Ad\j * Dfjt| + * (29) 

Only the strain increment S: A€ T + <A§>-<Ac> is responsible for T 
changing the size of the cylindrical volume, since the strain increment -Ac 
is used to force the cylindrical volume back to its original size after 
removal from the matrix. If the actual cylindrical grain with elasticity 
tensor J>f is now subjected to the strain increment J§: At + <Ae> - <Ac} it 
will aquire the same shape and size as the fictitious cylinder which has 
elastic properties identical to the matrix. The stress increment in the 
actual single crystal is then given by 

= 0t jk | { S k . TO tel * <**y>-< A6 n>} ' «<» 

The actual and fictitious cylinders now have the same final shape and size 
and if the stress increments in equations (29) and (30) are equal, the 
actual crystal cylinder can replace the transversely isotropic cylinder and 
still preserve continuity of displacements and tractions across the cylin- 
der-matrix interface. Equating (29) and (30) gives 

Aelj* - [( * < 31 ) 

Substitution of this result into equation (29) [or (30)] produces the stress 
increment in the cylinder due to the instantaneous application of the 
strain increment — <Ac> in the matrix in the form 

Affij « D-u (<**„> -<4e* i>) 

" Pjjti (3i) 

At the end of the creep response during the time interval At the 
stress increment in the single crystal cylindrical grain is given by equation 
(28). When the matrix is further elastically strained by the instantaneous 
application of the strain increment <As^ - the additional stress 
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increment in the single crystal grain is given by equation (31). Accor- 
dingly, the total stress increment in the single crystal due to the appli- 
cation of the strain increment in the time interval At is given 

by the Euler forward difference method as the sum of the stress incre- 
ments in equations (28) and (31), viz. 

x | (l*u- < ie t«>) -fcu' -< Ac t->) \ > (3J) 

where the first term in equation (31) is the definition of the matrix 
stress increment given in equation (17); 

Provided the stress level g is known in the single crystal at the 
beginning of the increment, the single crystal analysis in the preceding 
section furnishes the value of the inelastic strain increment Ac . Aver- 
aging bfi over the basal plane for crystals of different orientation by 
means of equation (18) then furnishes the value of . Assuming 

that the strain increment in the matrix is given, equation (17) furnishes 
the required stress increment in the matrix. The stress increment in the 
single crystal constrained by the transversely isotropic matrix is given 
by equation (32). All pertinent quantities can now be updated according 
to the Euler forward difference procedure and the process repeated for 
the next time increment. A procedure for evaluating the Eshelby tensor 
$ is given in the appendix. 


WORK IN PROGRESS 

The single crystal analysis has been coded into a FORTRAN sub- 
routine and is currently undergoing numerical test experiments. Coding 
of the directionally solidified analysis is due to commence shortly. Single 
crystal tubular specimens of PWA 1480 are being machined and will 
be tested at 1600°F under biaxial tension-torsion loading conditions. 
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APPENDIX 


The Eshelby tensor S is defined by the relation 

i5^riJ51, s tj(:-£')av( £ ') <» 

where <5ij(x - £ y ) is the elastic Green's function for the transversely 
isotropic matrix material. Although the Green's function for transver- 
sely isotropic materials is known [18], it is more convenient to work with 
the Fourier representation of the Green's function. The Fourier trans- 
form of . & ij ( X - r / ) is known [19,20] and it can be shown, by taking 
the Fourier transform of the defining equation, viz. 

and inverting the result, that 


6u(r-,')- — (ft d*K 


(3) 


In this equation the real part of the Fourier integral corresponds to the 
Green's function, and the Christoff el tensor M is defined by 


with 


«r o^s, 

£" K/Jk Iki = K/k 


(4) 


(5) 


being a unit wave vector in the direction of the Fourier wave vector £ . 

Introduction of this result into one of the integral terms in the deTini- 
tion of S gives, on reversing the order of the volume and wave vector 
integrations, a 


Lk *\) = M(r') 


(6) 


The volume integration extends over the cylindrical volume of the fictit- 
ious cylindrical grain and can be written as 

iOc^ + KiX-w+KjzO 


1 ' S/'-'.NCr') ■ du^ix, . 

Let x i w Then in cylindrical coordinates 

t* t* ( r" Uk, ( use-ntj^e) 

where a is the cylindrical radius. Since 

^ e 1 1 Ax^ s £ir , 

- 00 


(7) 


( 8 ) 


( 9 ) 
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where £( K* ) is the Dirac delta function, the integral takes the form 

I- a*SC«,)f f^"®****® 5 , dgd$. (10) 

Let g»<^/\J Kf+Kj * dgs * Then 

I- Mfcof -Jit ' 0.1 

•*o v° 

If we now set K \/{k*Zk} • toad' » then Ki/n/k^kJ sjig' , so that 


^ ** ^«0 $»0 


Since the integration extends over a whole circumference, it is imma- 
terial where the origin of 0 is placed. The integral may therefore be 
written as «£Cic> an 

I = i*«£ 5 J ifi 

«N«f vo Y-. 

. 2*SCK>) c*«W •» 

M tcf^t ) ^*v- a »wv 

° ^»0 

I* atofrtcj J, s (13) 

K, tK a v 

where T 0 and J, denote the usual Bessel functions. 

Equation (6) can now be written as 

. 1 -I&I _;K.r 

Wow *-K w \< # e ~~ so that 

k 1 

ws-i gs^f^s H-e? 

If k=3 or g=3, the Dirac delta function S( K$ ) gives zero values for 
the integral. Hence, the non-zero values of the Eshelby tensor S are 
given by k=l,2 and g s l,2. This arises because there are no compon- 
ents of the body force layer of surface tractions in the x s direction 
at the cylindrical interface between the grain and the matrix. Invoking 
the properties of the Dirac delta function gives 
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Tr .», \ 

Lk « 1 j 4Kl,lK,Mi j ft > 5 ‘’ t «* 0 ) tkC i 6 T^rJ, (oGfri) . 


2 TT 
|Iow put 




( 16 ) 


t,- T== * "* 6 » C»*- r=^= 
1 «|k,Sk{ J *,■>»<£ 


s si'*\6 , 


so that the Fourier wave vector components and K* correspond to the 
Z.i and sc* axes, respectively. Then in cylindrical coordinates, 

Lkfij — A 'J** 3-J^) 

or " a*© *«o K> 


l*« K*o 
at 


Uk(-;j »-J. f ev^j (t)c v 5 e ae jae‘‘ J|(«K)<K. 

1C*© 

Since S, is real, the real part of the preceding integral is 

X* f coos(Ki*)3\(*K)dK. 

«»o 

But, (cf. Gradstyn and Rizhik, p. 730, Eq. 6.671, No. 2) 

J- [a «S £ r 0£r£a. 

If 0*si li\'(r/a. ), then J* ae*sG/jo£7»- * e*se/[£CrVa]» ®sB/4*s6 * I . 
Thus 2c 

Lk rj " "in 1 ^ » 


(17) 


independent of position £ in the cylinder as expected from Eshelby's 
result. In this integral we have C,* cosS, £,=0> )” 1 

and k and g are restricted to the values 1 ana 2. The Eshelby 
tensor may now be written as 


SijW * DkjU ^ ^ ^ Mjj ^ 0*t j O^k d&J I' * 


(19) 


For a transversely isotropic material the Chistoffel tensor N has the 
component form 
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L 


D 



G 



The Eshelby tensor can now be evaluated by inverting the matrix in equa- 
tion (20) and integrating according to equation (19). Explicit results 

may be deduced, according to the calculations in Mura's book [21], for 
the resulting integrals. Checking of these integrals is currently in 

progress. 
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The twelve <110) slip direction vectors on the four octahedral {ill} 
planes. 



g. 2. Global and crystallographic axe- of the specimen. 




